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ABSTRACT 

Sh 2-188 is an example of strong interaction between a planetary nebula (PN) 
and the interstellar medium (ISM). It shows a single arc-like structure, consisting of 
several filaments, which is postulated to be the result of motion through the ISM. We 
present new Ha images from the Isaac Newton Telescope Photometric Ha Survey of 
the Northern Galactic Plane (IPHAS) which reveal structure behind the filamentary 
limb. A faint, thin arc is seen opposite the bright limb, in combination forming a 
closed ring. Behind the faint arc a long wide tail is detected, doubling the size of the 
nebula. The nebula extends 15 arcmin on the sky in total. We have developed a ‘triple¬ 
wind’ hydrodynamical model, comprising of the initial ‘slow’ asymptotic giant branch 
(AGB) wind and the later ‘fast’ stellar wind (the interacting stellar wind model), plus 
a third wind reflecting the motion through the ISM. Simulations at various velocities 
of the central star relative to the ISM indicate that a high velocity of 125 kms“^ 
is required to reproduce the observed structure. We find that the bright limb and 
the tail already formed during the AGB phase, prior to the formation of the PN. 
The closure of the ring arises from the slow-fast wind interaction. Most of the mass 
lost on the AGB has been swept downstream, providing a potential explanation of 
the missing mass problem in PNe. We report a proper motion for the central star 
of 30 ± 10 masyr“^ in the direction of the bright limb. Assuming the central star 
is moving at 125 ± 25 kms“^, the distance to the nebula is estimated to be 850('l42o 
pc, consistent with a spectroscopic distance to the star. Expansion velocities measured 
from spectroscopic data of the bright filaments are consistent with velocities measured 
from the simulation. Sh 2-188 is one of the largest PNe known, with an extent of 2.8 
pc. The model shows that this size was already set during the AGB phase. 

Key words: hydrodynamics - planetary nebulae: individual: Sh 2-188, S188, Simeiz 
22 - stars: AGB and post-AGB - ISM: structure - stars: mass-loss. 


1 INTRODUCTION 

The accepted theory of planetary nebula (P N) formation 
is the intera ctinK stellar winds model (ISW) JKwo 3 ll982t 
lBalicklll987ii where a fast wind (~ 10® kms ®) from the 
hot central star of a PN blows into the slow wind (~ 10 
kms“®) produced during the preceding asymptotic giant 
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branch (AGB) phase. The inner regions of the slow wind are 
compressed into a dense shell and ionised by the energetic 
UV radiation of the central star. The familiar ring-like ap¬ 
pearance of PNe is then observed. Structures in the nebula 
are normally attributed to asymmetries in the slow wind, 
related to physical properties of the central star, such as 
rotation or binarity. 

Observations of PNe have shown several cases where 
the outer shell shows the only departure from symmetry. 
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For those cases, the cause of the asymmetries has been pro¬ 
posed to be an interaction with the interstellar medium 
(I SM). Interac t ion of PNe with the ISM was first dis cussed 
by jGnrzadvanI 1^6^. An early theoretical study bv ISmithI 
assumed a thin shell a p proximation and the ‘snow¬ 
plough’ model of lOord lll95dl . ilsaacmamil lll97flll used the 
same approximation with higher velocities and ISM densi¬ 
ties. Both of these studies concluded similarly that a nebula 
fades away before any disruption of the nebular shell be¬ 
comes noticeable. 

In contrast, iBorkowski et all Jl99nll found that many 
PNe with large angular extent show signs of PN-ISM in¬ 
teraction, and that all nebulae containing central stars 
with a proper mot ion greater than 0.015 arcsecyr“^ do so. 
ISoker et al J (Il99lll hydrodynamically modelled the interac¬ 
tion. The PN shell is first compressed in the direction of 
motion and then in later stages this part of the shell is sig¬ 
nificantly decelerated with respect to the central star. Both 
conclude that the interaction with the ISM becomes dom¬ 
inant when the density of the nebular shell drops below a 
certain critical limit, of typically tih = 40cm“® for a PN in 
the Galactic Plane. 

IVillaver et alJ ll200,‘dl (hereafter referred to as VGM) 
pointed out the PN-ISM interaction had only been stud¬ 
ied by considering the relative movement when the nebu¬ 
lar shell had already formed. They performed 2D hydro- 
dynamic simulation s following the full AGB phase fo llowed 
by the PN phase iVassiliadis fc WoodI 119931 . Il99'dl . with 
a conservative relative velocity of the central star of 20 
kms“^ and conservative conditions of the surrounding ISM 
of riH = 0.1cm“®. VGM concluded that interaction pro¬ 
vides an adequate mechanism to explain the high rate of 
observed asymmetries in the external shells of PNe. Fur¬ 
ther, they conclude that stripping of mass downstream dur¬ 
ing the AGB phase provides a possible solution to the prob¬ 
lem of missing mass in PN whereby only a small fraction 
of the mass ejected during the AGB phase is inferred to be 
present during the post-AGB phase. Observational evidence 
for the effect of the ISM on AG B wind stuctures was found 
bv IZiilstra fc WeinbergeJ (|20^. 

The PN Sh 2-188 JSharDlesslll959ll is among the most 
extreme examples of ISM interaction. The nebula has a one¬ 
sided (semi-circular), filamentary appearance. It is a large 
nebu la, with a reported 340 arcsec diameter i Acker et alJ 
^93); the currrent paper shows it to be considerably larger. 
It is located in the Galactic plane at I = 128'’, b = —4°. 
New data are presented showing the faint back of the shell 
and an extended Ha tail. The unusual appearance suggests 
a high proper motion and makes Sh 2-188 an important 
test case for PN-ISM interaction at high velocity. We have 
developed a ‘triple-wind’ model using a initial slow AGB 
wind, a subsequent fast post-AGB (PN) wind, and adding 
a third wind reflecting the movement through the ISM into 
the ISW model. W e use a hydrodynamic scheme developed 
bv IWareind ||200^, to investigate whether this triple-wind 
model can reproduce the nebular shape of Sh 2-188 with¬ 
out requiring magnetic fields. We support the results of the 
model with a proper motion study of the candidate central 
star. 


Table 1. Estimates of the distance to, radius of and age of Sh2- 
188. 
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2 OBSERVATIONS 
2.1 Sh 2-188 

Sh 2-188 (PNG 128.0—04.1, also known as Simeiz 22) 
was photograph ed in 1965 by |Gaz^^^hain|_jW6^ Spec¬ 
troscopic work dParkei rrozinskav^VFs^ovI ll 97ll : 

I.IohnsorJ ll975L iKwi^iT revealed the object to have 

an extremely high [Nil] A6583 / Ha A6563 line-ratio. 
Further spectroscopic da ta on the filamentary structure 
JPosado fc Kwit,terill982lj indicated filament densities of a 
few hundred H cm“® and found abundances similar to Peim- 
bert Type I PN, where the large over-abundances of He and 
N result from CN O enrichment o f the stellar envelope prior 
to shell ejection iPeimberj|l98l |>. Mapping a t radio wave¬ 
lengths lllsrael fc Fel1illl97fi : ISalter et al.ll19S4^ found only a 
weak source which appeared to have a flat (thermal) spec¬ 
trum. This observational evidence, along with the relatively 
bright optical appearance and PN abundances, constituted 
a convincing argument to classify the object as a PN, rather 
than a supernova remnant where non-thermal radio emission 
wou ld be expected . 

Ik witter et alJ (ll988^ identified the candidate central 
star of the nebula with geometric, colour and apparent mag¬ 
nitude methods. The star is displaced from the geometric 
centre of the nebula. They estimate a percentage proba¬ 
bility of finding an unrelated blue star in the nebula to 
be 1-2 per cent at most. In t heir spectroscopic investiga¬ 
tion o f central stars of old PNe. lNapiwotzki fc Schoenbernerl 
tl995|) classified the c andid ate central star as hydrogen- 
rich. IBorkowski et al.l Jl99fll suggested that the filaments 
are most probably thin, sheetlike layers of shock-compressed 
nebular gas, parallel to the line of sight. The majority of the 
filaments are aligned along the crescent-shaped outer rim, 
but they noted a few appear to be along the nebular symme¬ 
try axis, possibly indicating stellar motion in that direction, 
as similar filaments do in Abell 35. 

Tabled shows various recent estimates of the distance 
and physical parameters of the nebula . Vario u s hist orical 
estimates can be found in table 5 ofJSaure^JjWQ^. We 
do not list statistical distances fe.g. iGahn et alJll992ll : the 
Shklovski met hod (co n stant PN mass) particularly is of too 
low accuracy. ISaureil (ll995^ calculated an extinction dis¬ 
tance based on published UBV photometry of stars in the 
angular vicinity used to construct a reddening distance re¬ 
lationship for the line of sight towards the nebula . Most 
recen tly, the NLTE modelling method of iNapiwotzkil lll999l . 
1200 has implied the nebula is ancient with the central 
star at the turning point before the UV flux switches off 
and the star enters the white dwarf cooling regime. The pa- 
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Figure 1. An image of Sh 2-188 created by a mosaic of images 
taken as part of the IPHAS survey in 2003 (see text). The faint 
central star is approximately 0.5 arcmin East of the brightest star 
in the centre of the nebula. North is up and East to the left. The 
image is 9.0 x 8.5 arcmin. 



Figure 2. Same as Fig. [Tlbut with scaling chosen to reveal faint 
structure. North is up and East to the left. The image shows a 
faint ring completing the bright semi-circular arc and tails stretch¬ 
ing away connecting in the north-west. The ring is centred approx- 
imately at RA = 1*' 30™ 25=, Dec = +58° 25' 30" (J2000). The 
image is 18.0 X 17.0 arcmin. 


rameters of Sh 2-188 are still subject to large uncertainties. 
Consistently though, this nebula is found to be an evolved 
object and as such is incl uded in the Atlas of Ancient PNe 
iTweedy_^J<witte3ll99 ^ . 

It weedy fc KwitteJ lll99^ note the existence of faint 
structure to the north-west of the bright nebula. They favour 
a simple explanation that the gas in this direction has dif¬ 
fused away. They infer a highly inhomogeneous ISM to the 
south-east producing the bright filaments but essentially ab¬ 
sent to the north-west. A more likely model is one of a high 
peculiar velocity through a uniform ISM, and this is explored 
below. 


2.2 Ha observations 

Fig .Qshows an image of Sh 2-188 created by a mosaic of new 
Ha A6563 observations taken as part of the Isaac Newton 
Telescope Photo metric Ha Survey of the Northern Galactic 
Plane (IPHAS) JPrew et alJl2nn,^ . The IPHAS survey en¬ 
tails wide-held imaging through Ha, r' and i'-band hlters. 
It uses the Wide-Field Camera (WFC) at the 2.5-m INT, 
located at La Palma. The goal of the survey is to image 
the entire northern Galactic plane between —5° and +5°, to 
study emission-line stars and nebulae. The WFC consists of 
four 4kx2k CCDs, in an L-shape conhguration. The pixel 
scale is 0.33 arcsec and the instantaneous held is about 0.3 
square degrees. Exposure times are 120 s for Ha and 10 s for 
i' and r' for the observations presented here. Each pointing 
is supplemented by a second one offset by 5 arcmin SW, to 
cover the gaps between the CCDs. The pointings together 
provide a tiled pattern that minimizes the gaps, but because 
of the L-shape CCD conhguration small gaps do appear. 

The area around Sh 2-188 was observed on 2003 Oct 13 
and 14. The seeing varied between 1 and 1.5 arcsec. The 
airmass was between 1.2 and 1.4. Data reduction was done 
using the pipeline procedures, which correct for the bias, 
hat-helding, non-linearity and ge ometric distortions . A de¬ 
tailed description can be found in lPrew et al.l (l2005h . Fields 
around the position of Sh 2-188 were retrieved and mosaiced 
together using the Virtual Observatory software package 
MONTAGE. The pixel scale was binned to 1 arcsec. A 40' 
held was mosaiced. The r' was subtracted from the Ha im¬ 
age to bring out the faint emission in this crowded held. 
MONTAGE tends to leave gradients across the full held, 
an artifact of the process of background matching of helds 
taken with different sky brightness. This linear gradient was 
subtracted. 

The hnal image in Fig. Q shows the bright portion of 
the nebula in the south-east, which describes a semi-circular 
arc with a diameter of 10 arcmin. Due to seeing variations 
between the r' and Ha frames, small core-halo residuals can 
be seen at the positions of bright stars. Stars can also leave 
positive or negative residuals due to the stellar absorption 
features in the Ha and/or r' hlter. The bright residual star 
just behind the brightest part of the shell is an M5 giant, 
where the r' filter covers a TiO absorption band. We stress 
that this is not an emission-line star, nor is it the central star 
of the PN. The central star of the PN is approximately 0.5 
arcmin East of this star. Negative residuals tend to identify 
stars with strong Ha absorption. The rectangular hole is one 
of the gaps left by the tiling pattern. 

IPHAS pr ovides both i mages and a source catalogue, 
as described in iDrew et alJ (l2005ll . The catalogue contains 
stars and small resolved objects but not large scale nebu¬ 
losities. Coordinates are derived by globally matching the 
stellar sources with all-sky astrometric catalogues. The in¬ 
ternal precision is within 0.1 arcsec, and agreement with 
external catalogues is 0.25 arcsec (e.g. USNO) to 0.1 arcsec 
(2MASS). We will use the IPHAS source catalogue below in 
the discussion on the proper motion of the central star. 

Structure first noted by iTweedv fc Kwitterl lll99(i) to 
the north-west of the nebula is now revealed by these ob¬ 
servations to be a faint circular closure of the bright nebula 
with connecting tails stretching behind. Fig. has scaling 
chosen to reveal this structure. From the front of the fila- 
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mentary structure to the connecting tails the nebula is 15 
arcmin long. The surface brightness of the tail is less than 
1 per cent of that of the bright filaments. The tail seems to 
be stretching away from the nebula in opposition to an in¬ 
ferred direction of motion which could explain the one-sided 
brightening as an interaction with the ISM. 

Fig. El shows an artificial combination of two three- 
colour images of the field. One is scaled to show detail in 
the bright filamentary structure and the other to show the 
faint structure. The r/ and the i' bands are coded as blue and 
green respectively, while Ha is coded as both red and yellow. 
The central star is indicated and is clearly displaced towards 
the bright arc from the geometric centre of the nebula. 


3 THE TRIPLE-WIND MODEL 

A hydrodynamical code has been used to simulate the inter¬ 
action with the ISM. The numerical scheme used to solve the 
hydrodynamic s equations uses the second-order Godunov 
scheme due to iFalla (Il99lh . The scheme is in 3D cartesian 
coordin ates, fully paral lel and includes the effect of radiative 
cooling iWareinsll200,‘jl . 

The numerical domain consists of a cubic grid with 200 
cells in each direction and uniform cell spacing. The central 
mass-losing star is placed at cell coordinates (50,100,100) 
and the simulation is performed in the frame of reference of 
the star. Mass loss is effected by means of setting the values 
of density, momentum and energy density in the cells in a 
spherical region of radius 5| cells centred on the star. The 
radius of the source volume has been chosen by an experi¬ 
mental process whereby it is a balance between producing 
the most spherical nebula when modelling the stationary 
ISW model with spherically symmetric winds (i.e. reduc¬ 
ing the pixelation of the cartesian grid) and not being large 
enough to interfere with the results of the simulation. The 
conditions within this region are reset at the beginning of 
every timestep to keep driving the wind. Movement through 
the ISM is parallel to the x-axis and material flows in at 
the [x = 1) boundary with a positive a:-velocity. All other 
boundaries have conditions allowing material to flow out of 
the domain freely. 

The hydrodynamic variables in the source volume are 
set according to the ISW model of the PN-ISM interaction. 
This follows the evolution of the AGB and post-AGB phases 
from the beginning of the AGB phase through the post- 
AGB or PN phase. In the triple-wind model, this evolution 
is combined with a constant movement through the ISM. 

Mass loss via a stellar wind has been modelled as a 
spherically symmetric constant mass-loss rate M with con¬ 
stant velocity v. In the ISM wind the density and velocity 
are constant. 

Typical estimates from the literature are used for the 
wind parameter values. The post-AGB fast wind parame¬ 
ters are: Mfw = 5 x 10“® M 0 yr“^, Ufw = 1000 kms“^ 
& Tfw = 5 X 10^ K. The AGB slow wind parameters are: 
Msw = 10“® MQyr”*^, Usw = 15 kms“^ & Taw = 10"* K. 
The switch between the AGB wind and the post-AGB wind 
is instantaneous and occurs after 10® years of AGB evolu¬ 
tion, which is suffiently long to allow a stable structure to 
develop. In view of the still considerable uncertainties on 
the detailed properties and evolution of these winds, more 


detailed temporal variations were not modelled. In reality, 
one may expect the AGB wind to show increasing mass-loss 
rates with time, whilst the post-AGB may increase in veloc¬ 
ity over time. However, the current assumptions are found to 
be sufficient to reproduce the basic structure of the nebula. 
The ISM itself is assumed to be homogeneous with charac¬ 
teristics of a warm neutral medium, the main consti tuent of 
the observed ISM: Tism = 2500 K & nn = 1 cm“® iBurtorJ 
^ 8 ^. The gas pressure in all three winds is calculated as¬ 
suming an ideal gas equation of state. 


4 RESULTS 

Various speeds of movement have been considered from 0 
hms”*^ to 175 hms”*^ in steps of 25 kms~^. T he results of 
all these models can be found in lWarein3 (l2005^ . We present 
here only the case of 125 kms“^ which we find best fits the 
general morphology of Sh 2-188. For comparison with obser¬ 
vations we have considered the locations of the brightened 
arc and its faint completing ring, the location of the can¬ 
didate central star (33 per cent of the geometric radius up¬ 
stream of the ring centre), and the location of the connecting 
tails downstream of the nebula. 

4.1 Post-AGB structure and evolution 

The left panel of Fig. El shows density on a slice through 
the computational domain at {y = 100 ) 1,000 years into the 
post-AGB (PN) phase. During the PN phase the wind de¬ 
creases in density and increases in velocity resulting in an 
adiabatic shock that forms the bright shell of the planetary 
nebula. Before this point, the nebula would have been inside 
the AGB wind bubble created by the AGB wind bow shock. 
At the depicted time, the nebula has expanded enough to 
be interacting with the AGB wind bow shock. Whilst the 
nebular shell has a constant temperature all the way round, 
the upstream part of the shell is now compressing the AGB 
wind bow shock and thus has a far higher density. Veloci¬ 
ties in this higher density region range up to 20 kms“^. The 
right panel of Fig. Elis a synthetic image obtained by squar¬ 
ing the density and collapsing the datacube along the line of 
sight. This is a simple approximation and does not account 
for temperature or optical depth. The nebula appears to be 
brightened upstream where the higher densities are located. 
The downstream portion of the nebula ring is faint yet clear 
behind the bright arc with the tails connecting together fur¬ 
ther downstream. 

Fig.Elshows the state of the PN-ISM interaction 2,000 
years into the PN phase. At this stage, the object is most 
reminiscent of Sh 2-188. The upstream arc is brightened, 
there is a faint completion of the nebular ring and the tails 
connect together downstream. Velocities in the upstream 
bright arc range from 20 to 70 kms“^. Importantly, the cen¬ 
tral star is now displaced 10 per cent of the geometric radius 
upstream of the geometric centre. The IPHAS observations 
of Sh 2-188 reveal the central star to be approximately 33 
per cent of the geometric radius upstream of the geomet¬ 
ric centre. The faint connection of the ring downstream is 
the part of the nebula responsible for this displacement. As 
the PN shell expands, it is slowed upstream by interaction 
with the AGB wind shock whereas downstream, the PN shell 
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Figure 3. An artificial combination of Fig.Qand Fig. [^showing the detail in the bright arc and the faint structure (with some blurring 
to allow combination of the images). The different IPHAS bands are colour-coded (see text). The Ha image includes the continuum. The 
faint central star is indicated on the image between the markers. 
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Figure 4. Shown on the left is a slice through the datacube showing the logarithm of density at a point 1,000 years into the PN phase. 
On the right is a representation of what the nebula might look like in emission. The images are 1 pc on a side. The scaling bar shown 
refers to the density slice and is presented in 10“® Mq pc“^. 









6 C. J. Wareing et al. 



IE4O4 2E404 3E404 4E404 5E404 6E+04 


Figure 5. Shown on the left is a slice through the datacube showing the logarithm of density at a point 2,000 years into the PN phase. 
On the right is a representation of what the nebula might look like in emission. The images are 1 pc on a side. The scaling bar shown 
refers to the density slice and is presented in 10~® Mq pc“^. 
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Figure 6. Shown on the left is a slice through the datacube showing the logarithm of density at a point 3,000 years into the PN phase. 
On the right is a representation of what the nebula might look like in emission. The images are 1 pc on a side. The scaling bar shown 
refers to the density slice and is presented in 10~® Mq pc“^. 


is progressing quicker through undisturbed AGB material, 
thus causing the geometric centre to move downstream. 

The structure of the PN-ISM interaction after 3,000 
years of the PN phase is shown in Fig. ini The nebula is now 
departing from circular symmetry and hence appearing less 
like Sh 2-188. The brightest regions are moving downstream 
to where the PN shell is still interacting with the shocked 
AGB material in the tails. The velocities in these regions 
range between 20 and 50 kms“^ as do the velocities in the 
shocked high density material at the head of the bow shock. 


The apparent displacement of the central star is approaching 
30 per cent of the geometric radius, similar to the star in 
the case of Sh 2-188. Eventually, after the UV flux of the 
central star has turned off and the star has entered the white 
dwarf cooling regime, the fast wind will cease. At this time 
the whole nebula is swept downstream, and the star begins 
to desert its nebula as it continues to move at high speed, 
leaving its tail behind. 
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Figure 7. Logarithmic density plots showing the development of the bow shock structure and tails during the AGB phase of evolution 
for a central star with a peculiar velocity of 125 kms“^. The images are Ipc on a side. Panel A (upper left) is after 12,500 years, Panel B 
(upper right) after 25,000 years. Panel C (lower left) after 37,500 years and Panel D (lower right) after 50,000 years. Density is presented 
in 10“® Mopc-3, 


4.2 AGB structure and evolution 

The formation and evolntion of the bow shock which so 
strongly influences the PN shaping is of interest. Fig.Qshows 
a slice of the computational domain through the position of 
the central star at various times in the AGB phase. The 
panels illustrate the development of the bow shock and tail 
structure. Note that the true structure can be thonght of 
as a rotation of the slices about the symmetry axis through 
the central star. Compression and deformation of the slow 
wind against the ISM rapidly forms a bow shock ahead of 
the star as shown in Panel A, 12,500 years after the start 
of the AGB phase. In Panel B, 25,000 years into the AGB 
phase, material stripped by ram pressure from the head of 
the bow shock is flowing downstream and beginning to form 
tails behind the nebula. Eventually, 37,500 years into the 
AGB phase, the flow collimates these tails downstream of 
the nebula and they appear to join. At this point in time 


the structure has reached a stable state with the position of 
the reverse shock ahead of the central star being predicted 
by a simple ram-pressure-balance argument. The temper¬ 
ature of the shocked ISM material at the head of the bow 
shock is approximately 200,000 K. This is in agreement with 
strong shock theory which predicts the temperature to be 
(3/16) (mvigyi/k), in this case a temperature of 205,500 K. 
Radiative cooling in the simulation has compressed this re¬ 
gion. The shocked AGB material is ten times more dense and 
at a temperature of a few thousand K. The bow shock ma¬ 
terial cools as it is ram-pressure-stripped downstream. The 
structure shown in the fourth panel of Fig. |7| is the state 
after 50,000 years of the AGB phase, with material clearly 
flowing back down the tails. This structure does not change 
for the next 50,000 years. The injected stellar wind balances 
the material moving downstream. 

An instability at the head of the bow shock causes ma- 
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terial to move down the tails in waves. This instability can 
be attributed to vortex shedding at the head of the bow 
shock, but the high flow velocity does not allow the vortices 
to form fully downstream. In our simulations of lower ve¬ 
locity central stars, this vortex shedding has been seen to 
destroy the smooth tail structure. This effect is not in the 
scope of this paper and will be discussed at length elsewhere. 

4.3 Proper motion of the central star 

Our simulation reproduces the observed features as de¬ 
scribed above, if the central star has a peculiar velocity of 
125 kms“^. If the star is within 1 kpc as predicted by vari¬ 
ous authors (see Table 0 and moving as we predict at 125 
kms“^, it should have an appreciable angular shift on the 
sky over a relatively short epoch. The heliocentric radial ve¬ 
locity is reported as Vr = —26 ± 10 kms“^. Corrections for 
the Solar motion and Galactic rotation approximately can¬ 
cel. Thus, a tangential velocity vt = 122 kms“^ with respect 
to the Local Standard of Rest (LSR) would be expected. 

With this in mind, we have performed a proper mo¬ 
tion analysis of the central star via astrometric methods. 
The observations used for this are detailed in Table Ini¬ 
tially we looked for a shift between the I PHAS observations 
discus sed earlier and the observations of iT weedy fc Kwitteil 
il99(tl taken in late 1992 at the 36 inch Burrell Schmidt 
telescope on Kitt Peak, using unbinned IPHAS images with 
a pixel size of | arcsec. Tweedy and Kwitter’s observations 
have a pixel size of 2.05 arcsec, which was insufficient to 
detect a shift in the position of the candidate central star. 

Following this we fou nd that the central star was listed 
in the USNO-Bl catalog llMonet et akllioo^ (star number 
1484-0053964). This catalog obtained data from scans of 
7,435 Schmidt plates taken for various sky surveys over the 
past 50 years. The candidate central star had been observed 
in 1976.9 to a high positional accuracy. We htted the po¬ 
sitions in the catalog to stars around the nebula up to a 
radius of 10 arcmin. The rms of the ht was 0.04 arcsec. We 
detected a signihcant shift of the candidate central star of 
0.8 arcsec. This shift equates to a proper motion of the cen¬ 
tral star of Sh2-188 of 26 ± 11 masyr“^ in right ascension 
and of —15 db 10 masyr“^ in declination. The 1992.8 epoch 
position encompasses the expected position of the central 
star at that time. 

This proper motion corresponds to a total displacement 
in the plane of the sky of 30 ± 10 masyr“^, with a position 
angle of 120° ± 20° East through North. The proper motion 
corrections due to Galactic rotation and due to the Solar 
motion approximately cancel in this direction. Assuming the 
central star is moving at 122±25 kms“^, the distance to the 
nebula is estimated to be 850jl42o P^- Th is distance r a nge is 
in go od a greement with the estimates of lNapiwotzkJ lll999l 
I 2 OOIII and lsaur'^ ^1995^ . The measured proper motion is in 
the direction expected from our simulation supporting the 
identffication of the central star of Sh 2-188. 

Although the agreement of our simulation with observa¬ 
tion is encouraging, this should not be overstated. We have 
investigated a large parameter range in velocity, but the 
range of possible density contrasts is very large and we do 
not claim the hnal result is unique. In fact, now the distance 
has been determined, we can also compare the physical size 
of the nebula with the simulation: this reveals that the ini¬ 


tial solution predicts a far smaller nebula than is observed. A 
solution to this discrepancy is discussed in the next section. 


5 THE SIZE OF THE NEBULA 

If, as the proper motion combined with the simulation pre¬ 
dicts, the central star is at 850 pc, there is a discrepancy over 
the physical size of the nebula between the observations and 
simulation. The observations show the nebula is 2.5 pc in 
diameter whereas in the simulation it is only around one 
tenth of this size. The large difference requires a change to 
the assumed values of the physical parameters. Since the 
proper-motion distance and other distance determinations 
agree, we keep the velocities unchanged and adjust the den¬ 
sities. This approach is also suggested by the fact that ob¬ 
served AGB (super) wind velocities vary by not much more 
than a factor of 2, whilst the mass-loss rates can vary by a 
factor of 100. A similar argument can be made for the ISM 
density versus the peculiar velocity. 

A simple rescaling of the simulation allows one to 
change the densities and radii, without changing the dimen¬ 
sionless model results. We wish to generalise the model hold¬ 
ing the velocities in the simulation constant and increase the 
size of the nebula by a factor of ten. The size of the nebula 
is dehned by the size of the AGB wind bubble which is in 
turn described by a simple ram pressure balance argument. 
To move the physical position of the ram-pressure balance 
point (reverse shock of the AGB wind bow shock) 10 times 
further upstream requires a factor of 100 increase from the 
ratio of slow-wind mass loss to ISM density, when keeping 
velocities constant. Increasing mass, length and time by a 
factor of ten decreases density by a factor of 100. Thus the 
simulation domain and nebula are ten times larger. The time 
scales also increase by a factor of 10: now 500,000 years on 
the AGB is required for a stable morphology to develop. This 
is in line with the 495,000 year AGB evolution predicted 
analytically bv IVassiliadis fc WoodI il993L [l9^ . The sur¬ 
rounding ISM has dropped in density from 1 cm~® to 0.01 
cm“®, a value consistent with the warm ionised medium 
(WIM) rather than warm neutral medium. The rescaling 
shows that we do not in any way predict the mass-loss rate 
in the slow wind, only the ratio between the mass-loss rate 
and ISM density. The snapshots of the nebula and AGB evo¬ 
lution presented in the results section become 10 times older 
and Sh 2-188 now appears most like the simulation after 20 
to 30,000 years of PN evolution. Thus, the distance, diam¬ 
eter and age are all in accordance with iNaniwotzkil lll999l . 
l2001^ and|s aure 3|T99i). 

Cooling times and cooling effects in the simulation are 
related to density and so decreasing the density by a fac¬ 
tor of 100 alters the cooling. We have rerun the simulation 
with this rescaling and found that the altered cooling pro¬ 
duces the same morphology as the initial simulation, thus 
the rescaling seems entirely reasonable. 


6 DISCUSSION 

The triple-wind model can reproduce the appearance of 
the nebula Sh 2-188. Accepting the model, the nebula has 
a peculiar motion with respect to the local ISM of 125 
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Table 2. A summary of observations of the central star of Sh 2-188 used in the estimation of proper motion. The A values give the 
measurement uncertainties. The IPHAS star has a magnitude of r' = 17.4 and r' — Ha = r' — i' = 0.0. 


epoch 

a. 

hh mm ss.s 

5 

dd ' " 

Aa 

n 

A5 

n 

Source 

1976.9 

01 30 33.113 

58 24 50.680 

0.128 

0.094 

USNO-Bl l^one^^t al. 2003') 

1992.8 

01 30 33.185 

58 24 49.540 

1.0 

1.0 

Tweedv k, Kwitter 119961 

2003.8 

01 30 33.203 

58 24 50.290 

0.165 

0.165 

IPHAS r-band IDrew et al. 20051 


kms“^. The peculiar motion is mostly in the plane of the 
sky with the radial velocity around -25 kms“^. The star 
is located close to the Galactic plane suggesting the star 
belongs to the disk population. The average transverse mo¬ 
tion of a P N-forming central star in the Galactic disk is 40- 
60 kms~ ^iBii:ine v fc Merrifleldlll99^ : IPgani fc SokeilllO^ 
ISkulian et alJll99% The tail of the distribution extends 

S I 130 kms~^ iSkulian*^ 'allll999t iFeast fc Whitelo^ 
. The velocity of the central star is therefore high, yet 
not unfeasibly high. It would be less extreme if the star were 
a thick disk star which have local av erage velocities of 50 — 75 
kms“^ jBinnev fc Merrifleldlll99lil . 

High peculiar velocities are more likely in older stellar 
populations. This is in conflict to the type-I classihcation of 
the nebula: the nitrogen enrichment indicates a more mas¬ 
sive progenitor. However, the high [N Ilj/Ha ratio may be 
caused by shock excitation, and does not necessarily imply 
an N overabundance. Further evidence of shock excitation 
is provided by the observation that the sulphur lines are 
stronger than H /l A4861 lines at all points where they have 
been measured liKwit,terlll97flll . 

The central star of Sh 2-188 belongs to the tail of the ve¬ 
locity distribution of (older) disk stars. Its rarity is borne out 
in the unique filamentary, one-sided appearance of the neb¬ 
ula, which is without parallel among known G alactic PNe. 
The simulations presented in IWarein j ^2005^ show that a 
velocity of less than 75 kms“^ does not result in connect¬ 
ing tails behind the nebula and between 50 and 100 kms“^ 
vortex shedding from the head of the bow shock seems to 
destroy the smooth tail structure observed in the nebula. 
A speed in the region of 125 kms“^ (using increments of 
25 kms“^) is thus the lowest speed from the simulations at 
which the morphology of Sh 2-188 is reproduced. 

The observed fragmentation of the bow s hock can be un¬ 
derstood in terms of the effect predicted by IPgani fc Sokerl 
il994Lir99ljI (hereafter referred to as PS). They postulated 
that Rayleigh-Taylor instabilities would fragment a bow 
shock in the direction of motion. Magnetic fields would sup¬ 
press certain modes of fragmentation and accentuate others, 
changing the appearance. Sh 2-188 is located in the Galactic 
plane and thus PS interpreted its fragmentation as being an 
effect of fluid and magnetic field instabilities. Importantly, 
PS conclude that fragmentation of a substantial part of the 
halo occurs only for a proper motion of greater than 100 
kms“^. This provides a lower limit to the proper motion of 
Sh 2-188 and supports the proper motion inferred from the 
simulation of 125 kms“^. There is no fragmentation in the 
simulation which could be attributable to a lack of resolu¬ 
tion, the lack of magnetic field or a lack of a gravitational 
field. Further, there is in reality temporal evolution of the 


slow wind, with thermal pulses at the end of the phase. It 
is possible that these pulses of high mass-loss could stratify 
and fragment the bow shock structure before the PN phase. 
As long as the lifetime of these stratifications allows them 
to survive late into the PN phase (which is questionable), 
the fast wind could reveal and further fragment them in a 
similar manner to the observed hlamentary structure. 

Spectroscopic data of the bright filaments taken by 
lR,osado fc Kwitt^ (ll982^ revealed velocities between 20 and 
70 kms“^. The velocities in the simulation show the same 
range in the high density regions at the head of the bow 
shock which in reality is thought to have fragmented form¬ 
ing the observed filamentary structure. This agreement be¬ 
tween the simulation and observations is very encouraging, 
although no more velocity data on this object are available 
to us for a more detailed comparison at this time. It should 
be noted that in our model these velocities are not related 
to the nebular expansion, but instead to the flow of shocked 
material from the head of the bow shock back down the tails. 

The faint emission completing the bright arc is particu¬ 
larly important. In the simulation this structure has a tran¬ 
sitory nature and moves quickly away from the central star: 
it traces the interaction between the fast post-AGB wind 
and the slow AGB wind (elsewhere in the nebula this front 
has merged with the ISM shock). One reason for the fast 
movement of this shock is the high ratio of post-AGB /AGB 
momentum which we use. 

The distance downstream of the ring from the central 
star compared to the distance upstream of the bright arc 
defines the shift of the geometric centre and provides a cal¬ 
ibrator for the simulation presented here. The rescaling of 
the simulation p redicts the age of the n ebula between 20 
and 25,000 years. iNaniwotzl^ (119991 12 OOill suggests that the 
star is 22,000 years into the post-AGB phase. 

The simulation is not entirely in agreement with the 
observations. The progressive displacement of a star from 
the geometric centre of a nebula is an indicator of evolu- 
tion and in the case of l arge nebulae, advanced evolution 
iTweedv fc Kwit,terlll99(t> . The central star of Sh 2-188 ap¬ 
pears to be approximately 33 per cent of the geometric ra¬ 
dius towards the bright arc from the geometric centre. The 
rescaled simulation indicates that after 20,000 years of PN 
evolution, the star should appear 10 per cent of the geomet¬ 
ric radius upstream of the geometric centre. 10,000 years 
later, the simulation shows the central star approaching 30 
per cent of the geometric radius upstream, in better agree¬ 
ment with the observations. However, by this time, the simu¬ 
lation predicts the bright regions of emission will have begun 
to shift downstream to where the PN shell is interacting with 
the tail. A better representation of the time-dependent stel- 
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lar winds m ay be needed to resolve this issue. iNaniwotzkil 
<199911200 ill locate the star near the start of the white dwarf 
cooling track. At this point, the wind is expected to cease. 
The forward ram shock will no longer be supported by in¬ 
ternal pressure and will retreat. The nebula will be blown 
downstream and disintegrate, and the white dwarf will ap¬ 
pear to move out of its nebula. This phase is not included in 
the simulation. However, it implies that the Sh 2-188 mor¬ 
phology has a limited life expectancy. 

Interaction with the ISM considerably alters the 
amount of mass within the circumstellar envelope during 
the AGB and post-AGB phases. The simulation reaches a 
point of stability during the life time of the AGB evolution 
llVassiliadis fc Woodlll99.'^ . Il99'dl : a longer AGB wind does 
not increase the mass of the shell. After 20,000 years of the 
PN phase, an investigation of the (rescaled) simulation re¬ 
veals that 0.6 Mq has been introduced into the spherical 
region of the nebula. This is made up of 0.5 Mq lost by the 
star and 0.1 Mq of ISM material entering the cross section 
of the nebula. An integration of the datacube over a roughly 
spherical region centred on the central star with the nebular 
radius estimates 0.236 Mq remaining within the PN bound¬ 
ary. Thus, 0.364 Mq of material has been swept downstream 
into the tail of the nebula. In fact, the observed tail of the 
nebula is far older than the nebula itself as it is made up 
of AGB and ISM material swept downstream. It is likely 
that most of the ISM material introduced into the stellar 
region has been swept downstream, but still a minimum of 
0.264 Mq of stellar material has been mixed with the ISM 
downstream providing an efficient method of returning stel¬ 
lar material directly to the ISM. 

Observations have revealed a missing mass problem in 
PNe whereby only a small fraction of the mass ejected dur¬ 
ing the AGB phase is inferred to be present during the post- 
AGB phase. Stellar evolution calculations predict that stars 
with initial masses in the range of 1-5 Mq will end as PN 
nuclei with masses around 0.6 Mq. Most of the mass is lost 
during the AGB phase and should be easily observable as 
ionised mass during the PN stage. However, observations 
of Galactic PNe have revealed on average only 0.2 Mq of 
ionised gas. Several PNe have been shown to have embedded 
molecular clumps which could contain much of this missing 
mass and these clumps have shown to survive ablation by 
stellar winds and th erefore could survive into the PN phase 
jPittard et alJl20^5^ . The stripping of mass downstream of 
a moving PNe as investigated here also provides an attrac¬ 
tive solution to this missing mass problem; the Dumbbell 
nebula appears to have clumps of material surrounding the 
bright nebu la inside what would a ppear to be an AGB wind 
bow shock llMeaburn et a DUq^ combining both of these 
solutions to the mi ssing mass problem. 

Observations jNapiwotzkil Il99'^ 1200 Jl have indicated 
that the central star of Sh2-188 has a mass of 0.6 Mq. 
These observations, in combination with the mass injected 
into our simulation and the addition of mass lost prior to 
the superwind phase (a typical estimate for low-mass stars 
being ^ 0.5 Mq) allows us to estimate a mass of ~ 1.5 Mq 
for the progenitor s tar. 

IVillaver et alJ ^2002^ discussed evolution of AGB stars 
and the interaction of their dust shells with the ISM. Typi¬ 
cally, AGB stars form large shells (up to 4 pc in diameter) 
around them. Evidence for such an AGB-ISM interaction 


front was found bv IZiilstra fc WeinbergeJ ||20o2), who found 
a M3 III AGB star surrounded by such a detached shell of 
4 pc diameter. The ISM influence is defined by the density 
and hence pressure in th e ISM, which drops off exponen¬ 
tially fr om riH = 2 cm~^ ^ISpitzer^ll978^ at a scale height of 
100 pc <Binnev fc Merrifiddr^O^ . The regular shells can 
be influenced strongly by the movement through the sur¬ 
rounding ISM. The AGB evolution of the central star of 
Sh 2-188 shows a structure extensively altered from these 
regular dust shells. Much of the AGB material has been 
swept downstream. 

At a distance of 850 pc, the nebula is 60 pc below the 
Galactic plane. We would expect the ISM density in the re¬ 
gion of Sh 2-188 to be around nn = 2 cm“® so it would 
seem our rescaling choice of nn = 0.01 cm“® is at odds with 
this. However, large regions of the Galactic plane in between 
spiral arms are dominated by su ch low-densit y gas, includ¬ 
ing the volume around the Sun l|Frischl ll995^ . In fact, the 
NE2001 model llGordes fc Lazioll20o3b used commonly in 
the Pulsar community to predict the galactic distribution 
of free electrons, predicts an electron density in the direc¬ 
tion of Sh 2-188 of rie = 0.013 cm“® at a distance of 850 
pc. This is in very good agreement with our prediction, but 
it should be stressed that this is only one choice of AGB 
mass-loss rate to ISM density ratio. A higher mass-loss rate 
would increase the ISM density, although this would have 
implications for the mass in the nebular region. The direc¬ 
tion of motion indicates the star is moving almost entirely in 
Galactic longitude and has travelled 64 pc during the AGB 
and post-AGB phases. Therefore the assumption of constant 
ISM properties during the AGB and post-AGB evolution of 
the central star may not be correct. 


7 CONCLUSIONS 

We have successfully reproduced the morphology and the 
available kinematic data of Sh 2-188 and understood its for¬ 
mation in terms of the ‘triple wind’ model. Eollowing the 
AGB evolution of the central star has been crucial in fitting 
the whole structure, in particular the tail behind the neb¬ 
ula which is comprised of purely AGB material. The triple¬ 
wind model of Sh 2-188 has predicted a velocity of the candi¬ 
date central star of 125 km s~^ . Velocities from spectrosc opic 
data on the bright filaments jRosado fc Kwitter^ll982^ are 
in agreement with velocities measured from the simulation. 
A proper motion study of the central star has shown it to be 
moving at 30.0± 10.0 masyr“^ in the direction of the head of 
the bright arc. The combination of these two measurements 
has resulted in estimates oi D — SSOI^jq pc, d ^ 2.5 pc and 
tpN = 22, 500 ± 2, 500 years. These estimate s are in agree¬ 
men t with the d i stance and age estimates of lSaureJ (^99^ 
and INaniwotzkil lll99fll . EoOlil . The prediction of ISM den¬ 
sity in the vicinty of Sh 2-188 is al s o in a greement with the 
NE2001 model of iGordes fc LazT^ (|20o3). The triple wind 
model explains the geometric displacement of the central 
star and indicates that the faint closure of the bright arc is 
a transitory structure which evolves downstream of the neb¬ 
ula. The PN-ISM interaction has caused « | of the mass 
expected in the region of the star to be swept downstream 
providing a solution to the missing mass phenomenon in PN 
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and a valuable way of mixing ISM and stellar material sev¬ 
eral pc downstream of the central star. 

The success of the triple-wind model to fit this extreme 
object gives confidence in our ability to fit objects with lower 
speeds. It is now clear that the outer halo structure of PNe 
contain the effects of ISM interaction and should not be 
modelled as stand-alone structures. Further, the ISM inter¬ 
action is an important method of mixing stellar material 
back into the ISM. The next generation of telescopes, par¬ 
ticularly ALMA, will be able to reveal cool dust structure in 
the universe and shed light on circumstellar AGB material. 
Further simulations considering temporal evolution of the 
stellar winds, magnetic fields and/or gravity may shed light 
on why we have not reproduced the observed fragmentation 
in Sh 2-188. 
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